Role of defect healing on the chirality of single-wall carbon nanotubes 
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Although significant efforts have been directed towards a selective single wall carbon nanotube 
synthesis, the resulting diameter and chirality distributions are still too broad and their control 
remains a challenge. Progress in this direction requires an understanding of the mechanisms leading 
to the chiral selectivity reported by some authors. Here, we focus on one possible such mechanism 
and investigate the healing processes of defective tubes, at the atomic scale. We use tight-binding 
Monte Carlo simulations to perform a statistical analysis of the healing of a number of defective 
tubes. We study the role of temperature as a primary factor to overcome the energy barriers involved 
by healing, as well as the role of the metal catalyst. Using both electron diffraction patterns and 
local characterizations, we show that the healing proceeds first along the tube axis, before spreading 
laterally, and observe the competition between two or more chiralities. The resulting picture is 
that no chirality seems to be favored by the healing mechanisms, implying that the reported chiral 
preference should result from other sources. 



PACS numbers: 61.46.Fg, 61.72.J-, 68.55.A-, 71.15.Nc 
I. INTRODUCTION 

Owing to their exceptional electronic properties, 
single- walled carbon nanotubes (SWNTs) stand among 
potential materials for future electronic and optoelec- 
tronic technologies. Although significant efforts have 
been put on synthesis, the nanotubes diameter and chi- 
rality distributions are still too broad and their control 
remains quite challenging. Sorting techniques ^ such 
as ultracentrifugation in a gradient density or mul- 
ticolumn gel chromatography ^ are presently able to se- 
lect nanotubes by diameters, chirality and also metallic 
or semiconducting character. However, such processes 
are complex and the tubes obtained are often defective 
and/or chemically spoiled. Such a structural control has 
also been reported to be achieved directly during the 
growtlPEl. It has been shown that the SWNT chirality 
distribution obtained on supported CoMo catalysts can 
be influenced by experimental conditions where a major- 
ity of near- armchair (6,5) can be produced^. Recently, 
a link between the composition of Ni^^Fei-a; nanocata- 
lysts and the resulting nanotube chirality has been ob- 
serveci^. These experimental results are not always well 
understood because the detailed microscopic mechanisms 
involved in the synthesis of SWNTs are still lacking. 

One does not know yet at which stage of the process, 
nucleation or growth, the diameter and chirality of the 
tube are fixed. One school of thought holds that the chi- 
rality is established during the nucleation stage. Static 
ab initio calculations have been performed to compare 
various forms of possible nuclei on the surface of metal 
and establish which cap or capped tube is energetically 
preferred^^ Reich et al. suggested that the nucle- 
ation of the cap fixes the chirality of an individual tube 
as a change in chirality is unlikely during the growth 



phas^I^. Yet, these calculations are performed on very 
crude atomic configurations of a cap sitting on a flat un- 
relaxed metal layer. In-situ transmission electron mi- 
croscopy (TEM) observations of the contact between the 
cap and the catalyst during the growth show a greater 
complexity suggesting that the chirality is established 
later during the growthPH^^. As an example, Hofmann 
et al. have observed in-situ the formation of a defected 
carbon nucleus from a metallic nanoparticlJ^, the de- 
velopment of a perfect tube being observed later. This 
implies that the defects observed at the initial stage have 
been annealed by some mechanisms, possibly involving 
the incorporation of additional C atoms into the net- 
work. Recent atomistic computer simulations by Neyts 
et al.^ using the so-called ReaxFF empirical potential, 
have confirmed this tendenc}'^^. By performing hybrid 
reactive molecular dynamics/force-biased Monte Carlo 
simulations, the authors found that the chirality of the 
growing SWNT cap can change, indicating that the ge- 
ometry of the tube is fixed later during the growth, in 
contradiction with the previous assumptions. 

Computer simulation is a unique tool to analyze the 
onset of the SWNTs chirality at the atomic scale. In- 
deed, the atomic resolution is not yet obtained in the in- 
situ TEM observation of the growtlP^HI^. Several mod- 
els have been proposed to study the nucleation-growth 
mechanisms of SWNTs^^ Nevertheless, wheth er the 
employed method is empiricaP^^^ or semi-empiricaJ^^^, 
all final configurations are plagued by a high concentra- 
tion of atomic-scale defects, i.e. non hexagonal rings, 
adatoms or vacancies. A specific attention must then 
be paid to the defect healing mechanisms, to establish 
whether they can be held responsible for the reported 
chiral selectivity. We have recently shown that defects 
in graphene are highly sensitive to the temperature and 
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discussed the key role played by metallic atoms in the 
reconstruction of a defected graphene sheet^^. 

In this paper, we investigate the healing processes of 
defective carbon SWNTs at atomic scale. We use a tight- 
binding model implemented in a Monte Carlo code to 
study the finite temperature evolution of defected tubes. 
Different lengths and diameters of nanotubes have been 
investigated at various temperatures. To understand the 
intrinsic character of the healing process we studied in- 
finite tubes, using periodic boundary conditions, and no 
cap. The onset of the tube chirality is investigated and 
discussed through both electron diffraction patterns and 
local analyzes. The approach proposed here enables one 
to try and identify the healing mechanisms that produce 
perfect tube structures and those possibly favoring a def- 
inite chirality. 



II. METHODOLOGY 

As in our previous works, a tight-binding potential is 
used to model the interaction between metal (here nickel) 
and carbon atoms. To keep the model as simple and fast 
as possible, the local densities of electronic states are cal- 
culated using the moment's method where only the s^p 
electrons of the carbon and d electrons of the metal are 
taken into account. This model, both simple and accu- 
rate, is then implemented in a Monte Carlo (MC) code 
using either a canonical or a grand canonical (GC) al- 
gorithm with fixed volume, temperature, number of Ni 
atoms, and C chemical potential /j^c- More details can be 
found in Ref.«^. Using our carefully parameterized order 
N tight-binding model, we have identified the tempera- 
ture and carbon chemical potential conditions for the nu- 
cleation of a tube cap on nickel particles. We showed that 
nucleation takes place after the outer Ni layer (s) are satu- 
rated with carborP^^. Recent technical improvement 
of the algorithm of our tight-binding Monte Carlo code 
made it significantly faster. The exploration of longer 
time scales and more extensive investigations are now 
possibl^^. 

We first generate a number of nanotubes with various 
realistic defects, by performing Grand Canonical Monte 
Carlo (GCMC) simulations where the active zone for 
inserting or removing carbon atoms is defined as a region 
of space around the surface of a cylinder (see Fig. [l^). 
Moreover, periodic boundary conditions have been 
applied along the z axis in order to avoid the closing of 
the tube. Indeed, due to the presence of dangling bonds, 
the spontaneous closing of the tube into a graphitic-like 
dome has been previously observed^^. The large number 
of defects, as shown in Fig. [T]3, results from the growth 
conditions: the carbon chemical potential chosen to 
obtain defected structures is unrealistically high (around 
-1.00 eV/at). Under such conditions, the growth is 
too fast and defects have a very high probability to be 
formed. In this way, we create a number of defected 
tubes with diameters and lengths ranging from 4.0 to 16 




FIG. 1: (a) Schematic representation of the cylinder area 
where GCMC are performed to form defected tubes, (b) Two 
examples of defected tubes considered in the present study 
and obtained after performing GCMC calculations at 2500 K 
and a value of /Jc about -1.00 eV/atom. 

A and 20 to 40 A , respectively. The initial configurations 
(153 to 756 atoms) contain in average 60 % of defects. 
The defects generally correspond to non hexagonal 
cycles and/or carbon atoms with a coordination number 
not equal to three. 

In a first step, we study the infiuence of the annealing 
temperature on the structure of a small tube with a di- 
ameter close to 4 A, as seen in Fig. It corresponds to 
the narrowest attainable tube synthesized, using AIPO4- 
5 as a templat d^^ l ^^ l The small tube, containing 180 C 
atoms, was subjected to temperatures ranging from 1000 
to 3000 K (see Fig. [2]3). In this case, no healing was ob- 
tained. On the contrary, the number of defects increases 
with temperature until linear chains are formed for tem- 
peratures around 3000 K. This is not surprising because 
such small tubes are only marginally stable, due to the 
strong curvature of the wall, and have been observed only 
in the presence of a template. 

III. HEALING MECHANISMS OF DEFECTIVE 
NANOTUBES 

A. Role of the temperature and chirality analysis 

We now focus on a larger defected nanotube whose 
diameter is between 6 and 8 A with periodic boundary 
conditions along the z axis (20 A). As in the previous 
case, this tube was subjected to temperatures ranging 
from 1000 to 3000 K in an attempt to heal the defects. 
The initial configuration contains 153 carbon atoms and 
60 % of the rings forming the tube walls are not hexagons. 
Figure |3] shows the final states of the MC runs at four 
temperatures (T= 1500, 2000, 2500 and 3000 K). At 1500 
and 2000 K, a significant amount of defects remains in 
the final configuration. Upon heating at 2500 K, the fi- 
nal configuration is an almost perfect tube where only 5 
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Monte Carlo Step (xlO'^) 

FIG. 2: (a) Defected tube with a diameter about 4 A. (b) 
Variation of the number of hexagonal rings as function of 
Monte Carlo steps and equilibrium configurations at 1000, 
2000 and 3000 K. 




FIG. 3: Equilibrium configurations at 1500, 2000, 2500 and 
3000 K. 



% of the defects remain. As seen in Fig. [4^, where the 
variation of the number of polygons during the MC sim- 
ulation is presented, 70 hexagonal rings are present after 
the healing process. At higher temperature (T= 3000 K), 
close to vaporization conditions, the tube is completely 
destroyed. The same conclusions hold for larger (diame- 
ters around 10 and 15 A) and longer (40 A) tubes studied. 
For systems without metal catalyst, 2500 K seems to be 
the optimal temperature, where the nanotube structure 
is relatively well healed. 

To go further, we have characterized the chirality of 
the healed tubes. The detection of armchair or zigzag 
tubes is straightforward, but other cases require a more 
sophisticated method: we chose to calculate their elec- 
tron diffraction patterns. This is a very powerful and 
reliable technique providing us with a direct determina- 
tion of the chirality. SWNTs are commonly defined by 



ir 

B 40 



(a) 


r\ 




6-Ring 


7-Ring 




5-Ring 













2000 4000 6000 8000 10000 
Monte Carlo Step (xio"^) 

(b) 



FIG. 4: (a) Variation of the number of polygons as function 
of Monte Carlo steps in a healed tube leading to a (6,4) at 
2500 K. (b) Diffraction patterns of tubes before (left) and 
after (right) healing process. 



the integer coordinates (n, m) of a so-called chiral vector 
Cn,m = na + mb in the 2D hexagonal lattice (a and b 
are unit vectors of the planar network). Given the di- 
ameter and the chiral angle of the tube, the chirality 
(n,m) can be determined. As an example, the compari- 
son of the electron diffraction patterns produced by the 
defected and healed tubes studied previously is depicted 
in Fig. [4^. All patterns have been obtained after perform- 
ing simulated annealing to relax the structures. This step 
is necessary to decrease as much as possible the elongated 
shape of the spots and then facilitate the interpretation 
of the images. In the present case, we have found two 
hexagonal networks rotated from each other by twice the 
chiral angle 10.5. Since its diameter is around 6.68 
A, we can assert that a (6,4) tube is observed after our 
healing process. 

To investigate how the tube's chirality builds up, we 
seek to identify the healing mechanism leading to almost 
perfect tubes. Our approach is to determine and fol- 
low the evolution of graphitic islands considered as seeds, 
with a local chirality based on the individual hexagons of 
the defected tubes. To distinguish a seed we have defined 
the following criterion. An area inside the defected tube 
is considered as a seed if at least four hexagons close to- 
gether are observed. This threshold has been established 
in an empirical manner from our statistical observations 
during our numerous simulations. The observed heal- 
ing process can be described as follows. Initially, the 



FIG. 6: Schematic representation of defect healing in the sim- 
ulations. 



FIG. 5: Equilibrium configurations at 2500 K of healed tubes 
(a) (5,5), (b)(6,3), and (c) (5,5) and (9,0). 



first seeds to be formed are destroyed by statistical fluc- 
tuations. When graphitic islands are sufficiently stable 
(Fig. |6|i), we have invariably observed their extension 
along the longitudinal axis of the tube leading to a rib- 
bon of hexagons (see Fig. [gJd). To confirm this mecha- 
nism, we have considered defected structures with a belt 
of hexagons all around the tube and investigate its be- 
havior at different temperatures. Instantaneously, most 
of the hexagonal rings forming the belt are destroyed and 
the process described previously can start again. Once 
the ribbons are formed (Fig. [g})), the growth of the em- 
bryo proceeds transversely, leading to a perfect tube as 
shown in Fig. To determine when the tube structure 
is established, a local chirality has to be defined, instead 
of the global one used above. Such analysis has been 
performed thanks to an approach proposed by Kim et 
al^. Using this method, for ideal SWNTs at K, the 
local chiral angles are exactly equivalent to the chiral an- 
gles of the tubes and the distribution histogram is single 
peaked, whereas at high temperatures it is broadened. In 
the case of defected SWNTs, the calculated local chiral 
angles show a multipeaked or broad distributions. If only 
one graphitic island with a defined chirality is present, 
the initial chirality is kept all along the process as shown 
in Fig. [7^. In the case where several seeds with different 
chiralities are identified, there is a competition between 
the formation of structures with one (Fig. ^}p) or many 
chiralities. The latter case is depicted in Fig. ^ where 
both (5,5) and (9,0) coexist until the end of the healing 
mechanism. 




FIG. 7: Variation of the chirality as a function of the Monte 
Carlo steps at 2500 K. (a) One initial chirality kept all along 
the simulation, (b) Competition between different seeds lead- 
ing to one chirality at the end of the healing process. 



B. Role of nickel 

We further investigate the correction and migration of 
defects in sp'^ structures by randomly adding Ni atoms 
on the walls of defective tube (Fig. [S^i). The aim of 
the present study is not to investigate the influence of a 



metallic nanoparticle in the growth process, which is very 
complex, but simply to understand the role played by 
metal atoms in the healing of defective carbon structures 
and also in the chiral selectivity. The system, containing 
153 C and 35 Ni atoms, was subjected to temperatures 
ranging from 1000 to 2500 K. As in the case of the pure 
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carbon structure, we notice that at 1000 K, the tube is far 
from being healed and a non-neghgible fraction of defects 
remains. Once again, large rings are healed at 1500 and 
2000 K, whereas pentagons and heptagons cancel out at 
higher temperatures. At 2500 K, the tube is completely 
healed as seen in Fig. [Sj This suggests that the healing of 
defects is thermally activated as observed for the isolated 
tube. By comparison with the case of tubes without cat- 
alyst, we deduce that Ni plays a crucial role in defect mi- 
gration. In fact, during the first steps of the simulation, 
we have observed Ni atoms preferring to cluster together 
instead of distributing over the entire SWNT(Fig. [sJd). 
A similar behavior has been noticed in our previous sim- 
ulation^ to explain the possibility of SWNT regrowth. 
As in the experiments where a vapor of metal atoms is 
deposited on the surface of SWNT seeds^^, we observed 
Ni atoms clustering at the lip of a rigid perfect tube to 
allow for regrowth. In the present instance, the nickel 
atoms tend to cluster and, when necessary, dissolve ex- 
cess carbon atoms to get the right number of atoms to 
form a better healed tube. Upon heating at 2500 K, the 
final configuration is an almost perfect tube where only 
2 % of the defects remain. The key role of Ni atoms 
in catalytically healing defects carbon structures has al- 
ready been observed with our TB model in the case of 
graphene^^ and by other grou ps in simulations on the 
growth of SWNT on Ni cluster?^. 

IV. FAVOURED CHIRALITY 

To investigate a possible chiral selectivity resulting 
from the healing mechanisms, we performed a study of 
the healing of twenty defected tubes taken as starting 
configurations: ten isolated tubes and ten tubes with 
nickel atoms dispersed on. According to their diameter 
range 6-8 A) and the periodic boundary conditions 
along the z axis (20 A), the possible chiralities were: one 
armchair, (5,5), two zigzags, (8,0) and (9,0), and three 
chiral tubes, (6,3), (6,4) and (8,2). Within our model, 
the cohesive energy differences between these tubes is 
negligible 1-2 meV/at), in agreement with ah initio 
calculation j^^ l ^^ l According to the calculated energies, 
no perfect tube is energetically favored with respect to 
others. After performing MC simulations at 2500 K, all 



defected tubes are healed. For 80 % of healed tubes, the 
exact determination of the chiralities through the calcu- 
lation of diffraction patterns is straightforward since few 
defects are still present. After the healing process, all 
possible chiralities were randomly observed. Some exam- 
ples of final equilibrium configurations are presented in 
Fig. |5j The same conclusions hold for larger and longer 
tubes. Our statistical study reveals that no chiral selec- 
tivity linked to the healing mechanism is to be expected. 
This is in agreement with TEM investigations where the 
atomic structure of several tubes was determined pre- 
cisely by combining a quantitative electron diffraction 
study and simulations based on kinematic theor}/^^. 
V. CONCLUSION 

To summarize, using tight-binding Monte Carlo sim- 
ulations, we have investigated the healing processes of 
SWNTs to determine how and when the chirality is es- 
tablished. We have shown that high temperatures condi- 
tions are capable of healing very defective tubes. Using 
this approach, we observed that strongly defected tubes 
are first repaired along their longitudinal axis, before the 
remaining defects are taken care of along the tube cir- 
cumference. Calculated electron diffraction patterns and 
local characterization of the chirality have been used to 
get a deeper insight into the onset of the SWNTs chiral- 
ity during their synthesis. A statistical analysis demon- 
strates the formation of tubes with uniform diameter but 
random chirality. The intrinsic defect healing mecha- 
nisms studied here, that do not include the presence of 
a catalyst nanoparticle at the lip, seem then unlikely to 
favor a specific tube chirality. Other phenomena occur- 
ring at the tube / catalyst nanoparticle interface, and in- 
fluencing the carbon atom incorporation energetics and 
kinetics, could play on the reported chiral selectivity and 
should be considered in future work. 
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